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Abstract 

Results are presented on a search for a light charged Higgs boson that can be pro- 
duced in the decay of the top quark t — > H + b and which, in turn, decays into t + v t . 
The analysed data correspond to an integrated luminosity of about 2 fb _1 recorded in 
proton-proton collisions at a/s = 7 TeV by the CMS experiment at the LHC The search 
is sensitive to the decays of the top quark pairs tt — > H ± W^bb and tt — > H ± H :p bb. 
Various final states have been studied separately, all requiring presence of a T lep- 
ton from H + decays, missing transverse energy, and multiple jets. Upper limits 
on the branching fraction £>(t — > H + b) in the range of 2-4% are established for 
charged Higgs boson masses between 80 and 160 GeV, under the assumption that 
B(H+ t+v t ) = 1. 
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1 Introduction 

The minimal supersymmetric extension of the standard model (MSSM) requires the introduc- 
tion of two Higgs doublets in order that the superpotential can contain appropriate terms for 
giving masses to both up and down type quarks HHSl- This leads to the prediction of five el- 
ementary Higgs particles: two CP-even (h,H), one CP-odd (A), and two charged (H ± ) states 
[9llT0|. The lower limit on the charged Higgs boson mass is 78.6 GeV, as determined by LEP 
experiments |TIT1 - IT4]| . If the mass of the charged Higgs boson is smaller than the difference be- 
tween the masses of the top and the bottom quarks, i.e. m H + < m t — m^, the top quark can 
decay via t — > H + b (charge conjugate processes are always implied throughout this paper). 
For values of tan /3 > 5, the charged Higgs boson preferentially decays to a t lepton and a neu- 
trino, H + — > t + v t , where tan j5 is defined as the ratio of the vacuum expectation values of the 
two Higgs boson doublets. In deriving the experimental limits we assume that the branching 
fraction £>(H + — > t + v t ) is equal to 1. 

The presence of the t — > H + b, H + — > t + v t decay modes alters the t lepton yield in the decay 
products of tt pairs compared to the standard model (SM). The upper limit on the branching 
fraction, B(t H+b) < 0.2, has been set by the CDF [15] and DO HU experiments at the 
Tevatron for m H + between 80 and 155 GeV, assuming £>(H + — > t + v t ) = 1. More recently, 
ATLAS experiment at the LHC has set the upper limit on the B(t — > H + b) between 5% and 1% 
for charged Higgs boson masses in the range 90-160 GeV [17J. 

The dominant process of production of top quarks at the Large Hadron Collider (LHC) is pp — > 
tt+X via gluon gluon fusion. The search for a charged Higgs boson is sensitive to the decays of 
the top quark pairs tf — > H+ W T bb and tf — >■ H+H+bb, where each charged Higgs boson decays 
into a t lepton and a neutrino. Throughout this paper, these two decay modes are referred to 
as WH and HH, respectively. 

Three different final states are studied, all requiring missing transverse energy and multiple 
jets. The T lepton decaying into hadrons and a neutrino is labeled T^. The first final state 
involves the production of Th and jets (labeled Th+jets), the second one is where Th is produced 
in association with an electron or a muon (labeled er^ or jiT^), and the third one is where an 
electron and a muon are produced (labeled ef/). Figure [ljshows representative diagrams for the 
Th+jets (left plot), e(^)ih (middle plot), and e}i (right plot) final states. We use a data sample 
recorded by the Compact Muon Solenoid (CMS) experiment until the end of August 2011 with 
an average number of interactions per crossing (pileup) of 5-6. The analyses correspond to an 
integrated luminosity ranging from 1.99 to 2.27 fb _1 depending on the final state. 

r + -> hadrons r + -> hadrons T+ — > £ + ^r'^ 




i r r 

Figure 1: Representative diagrams for the Th+jets (left), e(f/)Th (middle), and e}i (right) final 
states. 
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2 CMS detector, reconstruction, and simulation 



2 CMS detector, reconstruction, and simulation 

A detailed description can be found in Ref. ||T8| . The central feature of the CMS apparatus is a 
superconducting solenoid of 6 m internal diameter providing a field of 3.8 T. Within the field 
volume are a silicon pixel and strip tracker, a crystal electromagnetic calorimeter (ECAL), and a 
brass/ scintillator hadron calorimeter (HCAL). Muons are measured in gas-ionization detectors 
embedded in the steel return yoke of the magnet. Extensive forward calorimetry complements 
the coverage provided by the barrel and endcap detectors. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the x axis pointing to the centre of the LHC, the y axis pointing up (perpendicular to the LHC 
plane), and the z axis along the anticlockwise-beam direction. The polar angle 6 is measured 
from the positive z axis and the azimuthal angle (p is measured in the x-y plane. The preudora- 
pidity r\ is defined as — ln[tan(0/2)]. 

The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select the most interesting events in 
a fixed time interval of less than 4 ^s. The High Level Trigger (HLT) processor farm further 
decreases the event rate from around 100 kHz to around 300 Hz, before data storage. 

Muons are reconstructed [19| by performing a simultaneous global track fit to hits in the silicon 
tracker and the muon system. Electrons are reconstructed [20] from clusters of energy deposits 
in the electromagnetic calorimeter that are matched to hits in the silicon tracker. Jets, Th, and 
missing transverse energy (E™ 1SS ) are reconstructed using particles measured with the particle- 
flow algorithm Il2"l"| . The particle-flow algorithm reconstructs particles in each event, using the 
information from the tracker, the ECAL and HCAL calorimeters, and the muon system. Jets 
are reconstructed with the anti-fcx jet algorithm f22| with a distance parameter of R = 0.5. The 
value of E™ 1SS is defined as the magnitude of the vector sum of the transverse momenta of all 
reconstructed objects in the volume of the detector (leptons, photons, and hadrons). 

The b-tagging algorithm used in this analysis exploits as the discriminating variable the sig- 
nificance of the impact parameter of the track with the second highest significance [23|. The 
significance is defined as the ratio of the measured value of the impact parameter to the mea- 
surement uncertainty. The hadron-plus-strips (HPS) T identification algorithm [24J is used to 
reconstruct T leptons decaying hadronically. The HPS algorithm considers candidates with 
one or three charged pions and up to two neutral pions. The Th candidate isolation is based 
on a cone of AR = y 7 Acp 2 + At] 2 = 0.5 around the reconstructed Th-momentum direction. It is 
required that no charged hadrons with pj > pj ut and no photons with Ej > Ej ut be present 
within the isolation cone, other than the Xj, constituents. The typical values of pj ut and Ej ut are 
~1 GeV. 

Backgrounds tt, W+jets, Z+jets are generated with MadGraph 5 |25H26l interfaced with PYTHIA 
6.4.25 EZ). The diboson production processes WW, WZ, and ZZ are generated by PYTHIA. 
Single-top-quark production is generated with POWHEG |28|. The signal processes, tt — > H ± bH =F b 
and tt — > W ± bH :F b, are generated by PYTHIA. The TAUOLA 11291 package is used to simulate x 
decays in all cases. 

Generated events are processed through the full detector simulation based on Geant4 11301 131~1 , 
followed by a detailed trigger emulation and the CMS event reconstruction. Several minimum- 
bias events are superimposed upon the hard interactions to simulate pileup. The simulated 
events are weighted according to the measured distribution of the number of interaction ver- 
tices. The PYTHIA parameters for the underlying event were set according to the "Z2" tune, an 
update of the "Zl" tune described in Ref. fl32| . 
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The number of produced tt events is estimated from the SM prediction of the tt production 
cross section, 165lg(scale)^7(PDF) pb I|33ti36| . The theoretical prediction agrees with the cross 
section measured at the LHC H37ll5gl. 

3 Analysis of the r h +jets final state 

In the Th+jets analysis, events are selected by a trigger that requires the presence of a Th with 
transverse momentum pj > 35GeV and a large calorimetric EJj? lss (> 60GeV). The Th trigger 
selection includes the requirement on the leading-px track, pj > 20 GeV. The amount of data 
analyzed for this channel corresponds to an integrated luminosity of 2.27 ± 0.05 fb _1 . 

In this analysis, selected event are required to have one Th with p T h > 40 GeV within \tj\ < 2.1, 
and at least three other jets with pj > 30 GeV within \rj \ < 2.4 with at least one jet identified as 
originating from the hadronization of a b quark. 

In order to suppress the multijet background we use selection criteria on the missing trans- 
verse energy, E™ ss > 50 GeV, and on the angle between the E™ ss vector and the transverse mo- 
mentum of the Th, A^(p T h ,E T llss ) < 160°. This analysis selects T/, candidates with one charged 
hadron. The charged hadron is required to have Pj k > 20 GeV. In order to use non-overlapping 
data samples in the Th+jets analysis and the other analyses, events containing either an electron 
or a muon with p T > 15 GeV are rejected. The background events with W — > xv T decays are 
suppressed by a requirement on the variable R T = p trk / p Th , with R T > 0.7, which takes into 
account the different polarization of t leptons originating from H or W decays [39). Although 
the requirements on the transverse momenta of Th and the charged particle introduce a bias of 
R T requirement, it provides a background rejection factor of about two. 

In the Th+jets analysis the dominant reducible background arises from multijet events with 
large E™ ss and jets that mimic hadronic x decays or are misidentified as b-quark jets. 

The other background processes comprise electroweak (EWK) ones - W+jets, Z+jets, diboson 
(WW, ZZ, WZ) as well as SM tt and tW production. The W+jets and tt production processes 
dominate. These backgrounds can be divided in two parts: the first one labeled "EWK+tt x" 
consists of events where at least one x lepton in the final state is present with p T > 40 GeV, 
within \t] T \ < 2.1, and the second one labeled "EWK+tt no-T" consisting of events with no 
t leptons in the final state or with no x leptons satisfying the above-mentioned criteria. The 
"EWK+tt no-T" background events with no x leptons in the final state can pass the selection 
due to misidentification of a jet, an electron or a muon as a Th- 

The transverse mass, mj , can be reconstructed from the Th and E™ 1SS vectors, providing addi- 
tional discrimination between W and H decays. The shape and normalisation of the nij distri- 
butions of the multijet and "EWK+tt x" backgrounds are obtained from data. The nij distribu- 
tion of the multijet background is measured using the events which pass the signal selection 
described above, except for no requirements on x isolation and on an identified b quark jet. A 
small contamination from EWK+tt processes, evaluated using simulation, has been subtracted. 
The mj distributions are measured in bins of px of the x candidate (a x^ with no isolation cri- 
teria applied). The final nij distribution of the multijet background after full event selection 
is obtained by summing the nij distributions for each p T bin weighted with the efficiency that 
the t candidate passes the x isolation criteria and the R T selection. The efficiency is measured 
from data using events selected for the measurement of the mj distributions, but without ap- 
plying the requirements on Ei^ and A</>(p T h , E™ ISS ). The expected number of multijet events 
in a given bin i of the mj distribution is calculated as: 
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3 Analysis of the r h +jets final state 



^multijet _ ^multijet ^ p multi ) et £ T + R r (J) 
j 

where the index j runs over the bins of p T h ; £y +Rr is the efficiency of the x isolation and the 
R T selection in a given bin j, p™ ultl ^ et is the nij probability density function obtained from the 

shapes of the mj distributions, and ]V multl i et is the total number of the multijet background 
events. 

The expected number of events and the nij distribution of the "EWK+tt x" background are ob- 
tained using a control data sample defined with the same jet selection criteria of the Th+jets sam- 
ple, but requiring a muon instead of a T^. The reconstructed muons are then replaced by em- 
bedding in the events the reconstructed particles from simulated x lepton decays. The embed- 
ding method underestimates a small contribution from the Drell-Yan xx and WW— > xx + E™ ss 
processes, since a veto on the presence of a second lepton (e or }i) is used in the selection of the 
control sample. The residues of these backgrounds not counted with the embedding method 
have been estimated from the simulation. The "EWK+tt no-r" background has been estimated 
from the simulation. 

Figure [2] shows the event yield after each selection step starting from the requirement that at 
least three high-px jets are present. The expected event yield in the presence of the t — > H + b, 
H + — > x + v T decays is shown as the dashed line for m H + = 120 GeV and assuming £>(t — > 
H + b) = 0.05. The multijet background and the "EWK+tt x" background are shown as mea- 
sured from the data. The "EWK+tt no-r" background is shown as estimated from the simula- 
tion. 

The observed number of events after full event selection is listed in Table [TJ along with the 
expected number of events from the various backgrounds, and from the Higgs boson signal 
processes WH and HH at w H + = 120 GeV. The number of WH and HH events is obtained 
under the assumption that £>(t — > H + b) = 0.05. The systematic uncertainties listed in Table [j] 
will be discussed in Section [6] 

Table 1: Numbers of expected events in the Th+jets analysis for the backgrounds and the Higgs 
boson signal from HH and WH processes at m H + = 120 GeV, and the number of observed 
events after the final event selection. Unless stated differently, the expected background events 
are from simulation. 



Source 




HH + WH, m H+ = 120 GeV, B(t -+ H+b) = 0.05 


51 ±4 ± 8 


multijets (from data) 
EWK+tt r (from data) 

EWK+tt no-r 
residual Z/7* — > xx 
residual WW — > xv T xv T 


26 ± 2 ± 1 
78 ± 3 ± 11 
6.0 ± 3.0 ± 1.2 
7.0 ± 2.0 ± 2.1 
0.35 ± 0.23 ± 0.09 


Total expected background 


119 ±5 ±12 


Data 


130 



The nij distribution after all event selection criteria are applied is shown in Fig. [3] 
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tfs = 7TeV L = 2.3fb~ 1 CMS 
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x h +jets data 

- with hf-n^v 
I I multijets (from data) 
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Figure 2: The event yield after each selection step for the ih+jets analysis. The expected event 
yield in the presence of the t — > H + b, H + — > t + v t decays is shown as the dashed line for m H + = 
120 GeV and under that assumption that B(t -4 H+b) = 0.05. The multijet and the "EWK+tt 
t" backgrounds are measured from the data. The "EWK+tt no-r" background is shown as 
estimated from simulation. The bottom panel shows the ratio of data over background along 
with the total uncertainties. Statistical and systematic uncertainties are added in quadrature. 



4 Analysis of the er h and fir h final states 

The event selections used are the same as in the measurement of the top quark pair production 
cross section in dilepton final states containing T ||40l . 

In the eih analysis, the events are selected by a trigger that requires the presence of an electron, 
at least two jets with pj > 30 GeV and pj > 25 GeV, respectively, and a certain amount of H™ 1SS , 
where H^ 1SS is defined at the trigger level as the magnitude of the vector sum of pj of all jets in 
the event. As the peak instantaneous luminosity increased the requirements on the electron pj 
changed from 17 to 27 GeV and on H™ ss from 15 to 20 GeV. The amount of data analyzed for 
this channel corresponds to an integrated luminosity of 1.99 ± 0.05 fb _1 

In the f/Th analysis, the events are selected by a single-muon trigger with the threshold changing 
from 17 to 24 GeV during the data taking period. The amount of data analyzed for this channel 
corresponds to an integrated luminosity of 2.22 ± 0.05 fb _1 . 

The events are selected by requiring one isolated, high-px electron (muon) with pj > 35 (30) GeV 
and \rj\ < 2.5 (2.1). The event should have one Th with pj > 20 GeV within < 2.4, at least 
two jets with pj > 35 (30) GeV within \rj\ < 2.4, with at least one jet identified as originat- 
ing from the hadronization of a b quark, and E™ ss > 45 (40) GeV for the eTh (fTh) final state. 
The ih and the electron (muon) are required to have opposite electric charges. The isolation 
of each charged lepton candidate (e or }i) is measured by summing the transverse momenta of 
the reconstructed particles within a cone of radius AR = 0.3 around the lepton's direction. The 
contribution from the lepton itself is excluded. If the value of this sum divided by the lepton pj, 
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4 Analysis of the er h and pr h final states 



tfs = 7TeV L = 2.3fb~ 1 CMS 



120 GeV 



x h +jets data 

- with hf-n^ 
CD multijets (from data) 
CD EWK+tt i (from data) 

CD EWK+tt no-x (simul) 
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Figure 3: The transverse mass of Th and E™ ss after full event selection for the Th+jets analysis. 
The expected event yield in the presence of the t — > H + b, H + — > t + v decays is shown as the 
dashed line for m H + = 120 GeV and under the assumption that B(t — > H + b) = 0.05. The 
bottom panel shows the ratio of data over background along with the total uncertainties. The 
ratio is not shown for nij > 160 GeV, where the expected total number of the background 
events is 2.5 ± 0.3 while 5 events are observed. Statistical and systematic uncertainties are 
always added in quadrature. 



labeled J re i, is less than 0.1 (0.2) for electrons (muons), the lepton is considered to be isolated. 
The lepton is required to be separated from any selected jet by a distance AR > 0.3. Events 
with an additional electron (muon) with J re i < 0.2 and pj > 15 (10) GeV are rejected. 

The backgrounds in the eih and pTh final-state analyses arise from two sources, the first with 
misidentified Th, which is estimated from data, and the second with genuine Th, which is esti- 
mated from simulation. The misidentified Th background comes from events with one lepton (e 
or p), E™ 1SS , and three or more jets with at least one identified b quark jet (labelled "£+ > 3 jets" 
events), where one jet is misidentified as a t^. The dominant contribution to this background 
comes from W+jets, and from tt — > W + bW~b — > £vb qq'b (£ = e, p) events. The misidentified 
Th background is estimated by applying the probability that a jet mimics a Th to every jet in 
"1+ > 3 jets" events. The probability that a jet is misidentified as a Th is measured from data 
as a function of jet px and tj using W+jets and multijet events Ij24) . 

The backgrounds with genuine T leptons are Drell-Yan tt, single-top-quark production, di- 
bosons, and the SM tt events in which a t is produced from a W decay. The Z/7* — » ee, pp 
and tt -4 W + bW~b — > £ + vb£~vb events may also contain electrons or muons misidentified as 
Th- The event yields for these backgrounds are estimated from simulation. 

The data and the simulated event yield at various stages of event selection, described above, 
for the eTh (pTh) analysis are shown in Fig. [4] left (right). The backgrounds are normalized to 
the SM prediction obtained from the simulation. A good agreement is found between data and 
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the SM background. The expected event yield in the presence of t — > H + b, H + — > r + v T decays 
is shown as a dashed line for m H + = 120 GeV under the assumption that B(t — > H + b) = 0.05. 
The observed number of events after the full event selection is shown in Table [2]along with the 
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Figure 4: The event yields after each selection step for the eih (left) and fir^ (right) analyses. 
The backgrounds are estimated from simulation and normalized to the standard model predic- 
tion. The expected event yield in the presence of the t — > H + b, H + — > t + v t decays is shown as 
a dashed line for m H + = 120 GeV and under the assumption that £>(t — > H+b) = 0.05. The bot- 
tom panel shows the ratios of data over background with the total uncertainties. OS indicates 
the requirement to have opposite electric charges for a Th and a e or ]i. Statistical and systematic 
uncertainties are added in quadrature. 



expected numbers of events from the various backgrounds, and from the Higgs boson signal 
processes WH and HH for m H ± = 120 GeV. The misidentitied T background measured from 
the data is consistent with the expectation from simulation, 42 ± 4 (stat.) ± 8 (syst.) for the ei), 
analysis and 83 ± 7 (stat.) ± 12 (syst.) for the jir^ analysis. 

5 Analysis of the eji final state 

The event selections are the same as used in the measurement of the top quark pair production 
cross section in dilepton final states flU]. 

The e}i events are selected by a trigger requiring an electron with p| > 8 GeV and a muon with 
Pj > 17 GeV; or an electron with p| > 17 GeV and a muon with pj > 8 GeV. The amount of 
data analyzed for this channel corresponds to an integrated luminosity of 2.27 ± 0.05 fb -1 . 

In the eft analysis, the events are selected by requiring at least one isolated electron and at 
least one isolated muon (J re i < 0.15) in a cone of radius AR = 0.3 around the lepton with 
pj > 20 GeV within \rj\ < 2.5 (2.4) for electrons (muons). The event has to have at least two jets 
with pj > 30 GeV within \f]\ < 2.4. The leptons are required to be separated from any selected 
jet by a distance AR > 0.4. The invariant mass of electron-muon pair, m ef( , is required to exceed 
12 GeV. The electron and the muon are required to have opposite electric charges. 

The backgrounds considered in the e}i final-state analysis are the following: SM tt, Drell-Yan 
££ (£ = e, }i, t) production in association with jets (DY(££)), W+jets, single-top-quark produc- 
tion (dominated by tW) and diboson (WW, WZ, ZZ) production. Background yields are all 
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5 Analysis of the eji final state 



Table 2: Numbers of expected events in the er^ and fir^ analyses for the backgrounds and the 
Higgs boson signal from WH and HH processes at m H + = 120 GeV, and the number of ob- 
served events after the final event selection. Unless stated differently, the expected background 
events are from simulation. 



Source 


N|v h ± stat. ± syst. 


N|v h ± stat. ± syst. 


HH+HW, m H+ = 120 GeV, B(t -> H+b) = 0.05 


51 ± 3 ± 8 


89 ± 4 ± 13 


misidentified T (from data) 
tt -> WbWb -> ivb rvb 
tt -> WbWb — > £t/b £i/b 
Z/7* — > ee,}i}i 

Z/7* -> TT 

single top quark 
diboson 


54 ± 6 ± 8 
100 ± 3 ± 14 
9.0 ± 0.9 ± 1.8 

4.8 ± 1.8 ± 1.3 
17.0 ± 3.3 ± 3.0 

7.9 ± 0.4 ± 1.1 
1.3 ± 0.1 ± 0.2 


89 ± 9 ± 11 
162 ± 4 ± 23 
13.0 ± 1.2 ± 2.5 
0.7 ±0.7 ±0.7 
26.0 ± 4.3 ± 6.1 
13.5 ± 0.5 ± 1.9 
2.0 ± 0.2 ± 0.3 


Total expected background 


194 ± 8 ± 20 


306 ± 11 ± 32 


Data 


176 


288 



estimated from simulation. After the signal selection requirements are applied, 95% of the 
remaining background is due to SM tt decays. 

The data and simulated event yields at various stages of the event selection are shown in Fig. [5] 
The backgrounds are normalized to the standard model prediction obtained by simulation. A 
good agreement between the data and the standard model expectations is found. The expected 
event yield in the presence of t — > H + b, H + — > t + v t decays is shown as a dashed line for m H + = 
120 GeV under the assumption that <B(t —¥ H+b) = 0.05. It is smaller than the expectation from 
the SM alone (£>(t —¥ H + b) = 0) because the selection efficiency is smaller for H + — > t + v t — > 
£ + VfV T v T than for W + — > £ + Vf decay owing to the softer lepton pj spectrum. 

The numbers of expected events for the backgrounds and the Higgs boson signal processes 
from WH and HH modes at m H ± = 120 GeV, and the number of observed events after all 
selection requirements are summarized in Table [3] 

Table 3: Number of expected events in the eji analysis for the backgrounds, the Higgs boson 
signal from HH and WH processes at m H + = 120 GeV, and the number of observed events after 
all selection requirements. The expected background events are from simulation. 



Source 


Nev ± stat. ± syst. 


HH+WH, m H + = 120 GeV, B(t -> H+b) = 0.05 


125 ± 9 ± 13 


tt dileptons 


3423 ± 35 ± 405 


other tt 


23 ± 3 ± 3 


z/7* -> a 


192 ± 12 ± 19 


W+jets 


14 ± 6 ± 2 


single top quark 


166 ± 3 ± 18 


diboson 


48 ± 2 ± 5 


Total expected background 


3866 ± 38 ± 406 


Data 


3875 
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te = 7TeV L = 2.3fb 1 CMS 
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Figure 5: The event yield after each selection step for the e}i analysis. The backgrounds are 
from simulation and normalized to the standard model prediction. The expected event yield 
in the presence of the t — > H + b, H + — > t + v t decays is shown as a dashed line for m H + = 
120 GeV under the assumption that £>(t — > H + b) = 0.05. The bottom panel shows the ratios 
of data over background with the total uncertainties. The requirement for the e and ]i to have 
opposite electric charges is labelled as OS. Statistical and systematic uncertainties are added in 
quadrature. 



6 Systematic uncertainties 

The sources and the size of the systematic uncertainties are listed in Tables |4j |5j and [6] In all 
of the analyses the following effects are taken into account: 

• the uncertainty on the jet energy scale (JES), jet energy resolution (JER), and E™ 1SS 
scale. This uncertainty is estimated following the procedure outlined in Ref. [|42|: an 
uncertainty of 3% on the energy scale is included; 

• the theoretical uncertainties on the signal and background cross sections; 

• the uncertainty on pileup modelling due to the reweighting of simulated events ac- 
cording to the measured distribution of the number of vertices; 

• the uncertainty due to the limited number of events available in the simulated sam- 
ples (MC stat.); 

• an estimated 2.2% uncertainty in the integrated luminosity [43]. 

In addition, for the fully hadronic channel the following systematic uncertainties are taken into 
account: 

• the uncertainty on trigger efficiencies. The efficiency of the r part of the trigger is 
evaluated using Z — > tt events. It is used for the "EWK+tt x" background estimate. 
The data-to-simulation correction factor for the trigger on E™ 1SS is evaluated using tt 
events with an uncertainty estimated to be ~10%. The data-to-simulation correction 
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7 Evaluation of limits on B(t — ► H+b) 



factors for the efficiency of the trigger on Th and on E™ ss are used for the WH, HH 
signal and "EWK+tf no-r" background estimates; 

• the uncertainty on the estimate of the multijet background from data; 

• the uncertainty on the estimate of "EWK+tt x" background due to the uncertainty on 
the Th jet energy scale, the selection of muons in the control sample, the limited num- 
ber of events in the control sample, the contamination from multijet background, 
and the fraction of W — > x — > ]i events (/w^t^) in the control sample; 

• the uncertainty in the application of the lepton veto. It is estimated from the un- 
certainty in the lepton reconstruction, identification, and isolation efficiencies of 2% 
(1%) for electrons (muons), which is measured using Z — > ££ (£ = e, }i) events; 

In addition, for the analyses with Th in the final state (r^+jets, eTh, }i Th), the following systematic 
uncertainties are taken into account: 

• the uncertainty on the efficiency of x identification, estimated to be 6% [24]; 

• the uncertainty on the rate of misidentification of a jet as a Th or of a lepton as a Th, 
each estimated to be 15% 11241 : 

• the uncertainty on the efficiency of b tagging, 5.4% II23II : 

• the uncertainty on the rate of misidentification of a jet as a b quark, 10% Il23l : 

In the eTh and jix^ analyses the uncertainty in the estimation of the misidentified x background 
has two sources: the limited number of events for the measurement of the X misidentification 
rate and the difference in the X misidentification rates for jets originating from a quark with 
respect to jets originating from a gluon. 

Finally the uncertainty on the reconstruction, identification, and isolation efficiency of an elec- 
tron or a muon is taken into account in the eTh, and efi analyses. It is estimated to be 
-2-3%. 

The full sets of systematic uncertainties are used as input to the exclusion limit calculation. 

In the T+jets analysis the mj distribution shown in Fig. [3] is used in a binned maximum- 
likelihood fit in order to extract a possible signal. Other channels use event counting only for 
setting the limits. The uncertainties on the shapes for the multijet and "EWK+tt x" backgrounds 
derived from data are evaluated taking account of the corresponding uncertainty in every bin 
of the mj distribution. In addition, the mj shape uncertainty for the "EWK+tt x" background, 
related to the Th energy scale uncertainty, is taken into account in the fit. For the signal and 
the small "EWK+tt no-T" background the mj shape uncertainty in the JES+JER+E™ 1SS scale is 
evaluated from simulation. 

7 Evaluation of limits on B(t — > H+b) 

The expected number of tt events, after final event selection, is shown in Fig.|6]for the fix^ (left) 
and e}i (right) analyses as a function of the branching fraction B(t — > H + b) for m H + = 120 GeV. 
Expectations are shown separately for contributions from WH, HH, and tt — > WbWb (WW) 
processes. In the eTh, F T h/ an d fully hadronic analyses the total tt event yield ^^ SSM ) from 
WW, WH, and HH processes is larger than the yield from the standard model tt — > WbWb 
process (N? M ). This is due to the fact that the branching fraction for the Higgs boson decay into 
xv T is larger than the corresponding branching fraction for W boson decay. For the eji analysis 
the total tt event yield is smaller than that expected from the standard model. 
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Table 4: The systematic uncertainties on event yields (in percent) for the ih+jets analysis for 
background processes and for the Higgs boson signal processes WH and HH in the range of 
m H + = 80-160 GeV. The range of errors for the signal processes is given for the Higgs boson 
mass range of 80-160 GeV. 
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WH 


multi- 
jets 


EWK+tt t 


EWK+tt no-T 
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tt 
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W+jets 
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26 


23 
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cross section 


+7.0 
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+7.0 
-10.0 






4.0 


4.0 


+7.0 
-10.0 


8.0 


5.0 


pileup modeling 


0.3-4.2 
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7.6 


3.9 


7.1 


15 


10 


MC stat 


6.2-11 
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28 
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luminosity 
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11 


12 


11 


12 


11 
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15 
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1.6 




jet— > b misident. 










2.0 


2.6 






4.8 



Table 5: The systematic uncertainties on event yields (in percent) for the fiT^ analysis for the 
background processes and for the Higgs boson signal processes WH and HH for m H + = 
120 GeV. 
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tt« 


misident. r 
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DY( W ) 


DY(tt) 
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4.0 
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9.0 
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15.0 
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Table 6: The systematic uncertainties on event yields (in percent) for the eji analysis for the 
background processes and for the Higgs boson signal processes WH and HH at m H + = 
120 GeV. 





HH 


WH 


tt 


DY(U) 


W+jets 


Single top 


diboson 


JES+JER+E™ 88 


2.1 


2.0 


2.0 


6.0 


10.8 


4.0 


6.5 


cross section 


+7 
-10 


4.3 


5.0 


7.4 


4.0 


pileup modeling 


4.5 


4.5 


5.0 


5.5 


4.0 


5.5 


5.5 


MC stat 


5.3 


7.9 


1.0 


6.5 


42.9 


1.9 


4.3 


luminosity 


2.2 


dilepton selection 


2.5 




0.8 1 0.2 0.4 0.6 0.8 1 

e(t^H + b) 6(t^H + b) 



Figure 6: The expected number of tt events after event selection for the \it^ (left) and e}i (right) 
final states as a function of the branching fraction £>(t — > H + b) for m H + = 120 GeV. Expecta- 
tions are shown separately for the WH, HH, and WW contributions. 
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Assuming that any excess or deficit of events in data, when compared with the expected back- 
ground contribution, is due to the t — > H + b, H + — > r + v T decays, the value of x = B(t — > H + b) 
for each individual analysis can be related to the difference AN between the observed number 
of events and the predicted background contribution through the following equation: 



AN = N 



MSSM 
tt 



N. 



SM 



2x(l - x)N WH + * 2 N HH + [(1 - x) 2 - 1]N. 



SM 



(2) 



In this equation Nwh is estimated from simulation forcing the first top quark to decay to H ± b 
and the second to W T b, and Nhh forcing both top quarks to decay to H ± b. In the eih, MJh/ anc 
e}i analyses, N? M is evaluated from simulation, as given by the tt background in Table |2] and [3 
In the Th+jets analysis, most of the tf —¥ WbWb yield is derived directly from data, so it does 
not contribute to AN whatever the value of x. In other words if an H + SUSY signal is present in 
the data, affecting the tf — > WbWb rate, it also affects the data driven background estimate for 
this rate and therefore this contribution disappears in the difference data — background. In this 
case N? M contains only the small tt contribution included in the "EWK+tt no-r" background 
in TablejlJ which is derived from simulation: N? M = 2.1 ± 0.6 (stat.) ± 0.5 (syst). 

The CLs method ||44l [451 is used to obtain an upper limit, at 95% confidence level (CL), on 
x = £>(t — > H + b) using Eq. [2] for each final-state analysis and for their combination. The 
background and signal uncertainties described in Section [6] are modeled with a log-normal 
probability distribution function and their correlations are taken into account. In the T+jets 
analysis the mj distribution shown in Fig. [3] is used in a binned maximum-likelihood fit in 
order to extract a possible signal. For the eih, }IT^, and eji final states only event counting is 
used to obtain the upper limits. 

The upper limit on B(t — > H + b) as a function of m H + is shown in Fig. |7]for the fully hadronic 
and eTh final states and in Fig. [H] for the jiT^ and eji final states. The combined upper limit 
has been obtained using the procedure described in [46]. Figure [9] (left) shows the upper limit 
obtained from the combination of all final states. 



i/s = 7TeV L = 2.3fb~ n CMS 



^0.14 

T 

O 0.1 



- t->H*b, H*^tv 
x h +jets final state 
S(H*^xv) = 1 



■ Observed 

■ Expected median +1 a 

■ Expected median +2 a 



-O 

+ 

X 

T 

H — ' 

CQ 
O 



0.4 



0.35 



0.3 



\fs = 7TeV L=2.0fb" 1 CMS 



'0.25 




t^H*b, HWcv 
ex h final state 
B(H*^tv) = 1 



Observed 
' Expected median ± 
. Expected median ± 



130 140 150 160 

m H+ (GeV) 



130 140 150 160 

m H+ (GeV) 



Figure 7: Upper limit on B(t — > H + b) as a function of m H + for the fully hadronic (left) and the 
eTh (right) final states. The ±1(7 and ±2c" bands around the expected limit are also shown. 



Table [7] gives the values of the median, ±.\a, and ±2<7 expected and the observed 95% CL 
upper limit for £>(t — > H + b) as a function of ?«h+ for the combination of the fully hadronic, 
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0.4 



i/s = 7TeV L = 2.2fb" 1 CMS 



0.35 



+ X 

T 

H — ' 

CQ 0.3 
o 



1 1 1 1 1 1 1 1 
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Figure 8: Upper limit on <B(t — > H + b) as a function of m H + for the fiT^ (left) and eji (right) final 
states. The ±1(7 and ±2c bands around the expected limit are also shown. 



eih, jiTh, and e}i final states. The systematic uncertainties for the ei^, jir^, and eji analyses 
are larger than the statistical uncertainties. Figure [9] (right) shows the exclusion region in the 

Table 7: The expected range and observed 95% CL upper limit for B(t — > H + b) as a function of 
m H + for the combination of the fully hadronic, eih, pr^, and eji final states. 



95% CL upper limit onB(t^ H+b) 


m H+ 
(GeV) 


Expected limit 


Observed 
limit 


-2a 


-la 


median 


+la 


+2a 


80 


0.018 


0.022 


0.029 


0.040 


0.054 


0.041 


100 


0.014 


0.018 


0.024 


0.032 


0.043 


0.035 


120 


0.013 


0.015 


0.020 


0.027 


0.040 


0.028 


140 


0.009 


0.011 


0.014 


0.021 


0.030 


0.022 


150 


0.008 


0.010 


0.013 


0.019 


0.027 


0.023 


160 


0.008 


0.009 


0.011 


0.016 


0.023 


0.019 



MSSM m H +-tan f> parameter space obtained from the combined analysis for the MSSM m^ ax 
scenario EZI: M SU sy = ITeV, ji = +200 GeV, M 2 = 200 GeV, trig = 0.8M S tjsy, X t = 2M SU sy, 
and Ab = A t . Here, Msusy denotes the common soft-SUSY-breaking squark mass of the third 
generation; X t = (A t — ]il tan/S) is the stop mixing parameter; A t and Ay, are the stop and 
sbottom trilinear couplings, respectively; ji the Higgsino mass parameter; M g the gluino mass; 
and M-2 is the SU(2)-gaugino mass parameter. The value of M\ is fixed via the unification 
relation Mi = (5/3)M2sin0 w / cos# W - 

The t — > H + b branching fraction is calculated with the FeynHiggs program 11481 . The exclusion 
contours corresponding to the ±1(7" theoretical error on B(t — > H + b) due to missing one-loop 
EW corrections (5%), missing two-loop QCD corrections (2%) and Aj, induced uncertainties 
(the Afc term accumulates the SUSY-QCD corrections) [36 J are also shown in Fig. [9] (right). 

The upper limit on the the branching fraction B(t — > H + b) and the exclusion region in the 
MSSM m H +-tan /3 parameter space obtained from the combined analysis are comparable with 
the results from the ATLAS experiment fTT) . 
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m H+ (GeV) m H+ (GeV) 



Figure 9: Left: the upper limit on6(t — > H + b) as a function of m H + obtained from the com- 
bination of the all final states. Right: the exclusion region in the MSSM Mn+-tan j5 parameter 
space obtained from the combined analysis for the MSSM mj^ ax scenario Ii47ll . The ±lcr and 
±2cr bands around the expected limit are also shown. 

8 Summary 

A search has been performed for a light charged Higgs boson produced in top quark decays 
t — > H + b and which in turn decays into t + v t . The data sample used in the analysis corresponds 
to an integrated luminosity of about 2 fb _1 . The fully hadronic, eih, jit^, and eji final states have 
been used in the analysis. The results from these analyses have been combined to extract limits 
on t -> H + b branching fraction. Upper limits on the branching fraction H(t — > H + b) in the 
range of 2-4% are established for charged Higgs boson masses between 80 and 160 GeV, under 
the assumption that £>(H + — > r + v T ) = 1. 
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